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ABSTRACT 


The  accuracy  of  the  lower  confidence  limit  procedure  in 
NAVwEPS  OD  2930^  is  analyzed  wnen  the  failure  distributions 
of  the  components  in  the  series  system  are  either  gamma  or 
truncated  normal.   Several  representations  of  the  accuracy 
are  supplied. 
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CHAPTER  I 
INTRODUCTION 

The  statistical  model  developed  and  stated  in  NAV^EPS 
OD  29304  "Guide  Manual  For  Reliability  Measurement  program" 
was  developed  for  the  Department  of  the  Navy,  Special 
Projects  Office.   It  contains  very  versatile  procedures  for 
obtaining  both  system  reliability  point  estimates  and  system 
reliability  lower  confidence  limits  using  a  variety  of  test 
data  on  component  parts,  subassemblies,  assemblies,  and 
subsystems. 

This  document  makes  the  following  assumptions  on  the 
components  of  the  series  system: 

(1)  Constant  failure  rate 

(2)  Additivity  of  stress  effects 

(3)  Independence  of  part  failure 
{k)  Failure  rate  constancy 

The  purpose  of  this  thesis  is  to  check  the  accuracy  of 
this  statistical  model  without  the  assumption  of  a  constant 
failure  rate.   The  model  will  be  checked  for  a  system  com- 
posed of  four  components  in  series,  using  various  sample 
sizes,  mission  times,  and  combinations  of  the  gamma  and 
truncated  normal  failure  distributions.   The  statistical 
model  will  be  analyzed  by  computer  simulation  using  a  Monte 
Carlo  method  to  generate  the  gamma  and  truncated  normal  test 
times.   The  parameters  of  the  distributions  were  selected  so 


that  the  failure  rate  curves  would  be  representative  of 
those  curves  used  by  contractors  or  that  the  system  relia- 
bility was  a  preselected  value. 
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CHAPTER  II 

THE  STATISTICAL  MODEL 

Section  three  of  reference  1  contains  a  detailed 
explanation  of  the  statistical  model  which  is  being  analyzed. 
The  following  definitions  are  applicable  parts  of  the  model 
used  in  the  computer  simulation  for  the  lower  confidence 
limit  of  a  series  system  with  four  components.   For  a  given 
sample  size,  N,  ,  i  =  1,2,3,^1-  and  a  predetermined  test  time, 
Toi ,  i  =  1,2,3^  a  best  estimate  of  the  failure  rate  for 
component  i  is  given  by: 

f,   2N, 

1  =  „  1  ?„T  \  .  (1) 

-i    Sj   2^  +  1 

l 
Where  Si  =  £  T*  1 ,  is  the  sum  of  all  test  times 

accumulated  on  the  N  components  of  type  i,  and  f.  is  the 

number  of  failure  times  which  are  less  than  T  . .   From  these 

oi 

best  estimates  of  failure  rates,  the  system  reliability 
estimates  are  given  by: 

Rs  =  e  (2) 


where  X  =   2  k>  (3) 

s    1  =  1  1 


The  model  then  uses  the  following  formulae  in  order  to 
compute  the  lower  confidence  limit  on  system  reliability: 
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A    2\    +  (K,3)2+  V  *Ua(K|3)2C+  (K^C2 

Au  =  _S 1 (4a) 


k 
if   2  f  j_  >  0 
i  =  1 


J  =IMI!   ^   JL  (,b) 

u     n    i  =  1  S± 


k 
if    £  f<  =  0 

1  =  1 


where  n  is  the  number  of  tests. 

a 

Here  \  is  the  upper  confidence  limit  of  the  failure 

rate  for  a  system  with  k  units  in  series,  and  K  is  a  per- 
centage point  of  the  normal  distribution.  3  is  a  correction 

factor  for  small  values  of  \.    and  c  =   2  — /  A._.   All  times 

1         1  =  1  si   S 


are  in  mission  units. 


A 


Reference  1  then  defines  R  T  /-  \  to  be  the  lower  con- 
fidence  limit  on  system  reliability: 

-K  ' 

Es,L(a)  =  e  (5) 

It  must  be  noted  that  since  equation  (2)  is  only  a  best 
estimate  for  the  failure  rate,  the  statistical  model  ex- 
plained serves  only  as  an  approximation  for  system  relia- 
bility. 
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CHAPTER  III 


SIMULATION  PROCEDURE 


The  system  to  be  simulated  will  consist  of  four  com- 
ponents in  a  logical  series,  with  a  system  reliability,  R  , 

s 

which  may  be  expressed  as: 


Hs  =   I  I   Hi 

1  =  1 


wnere  R  =  P(T>1)  is  the  true  reliability  of  the  i  Xi  com- 
ponent at  mission  time  one.  It  is  desired  to  find  a  lower 
100(1  -cl)%   confidence  limit  for  R  ,  namely  the  statistic 

R  T(    ^.      If  S    .  .  is  in  fact  the  exact  100(1-  a )%   lower 
s,L(a)       s,L(a) 

confidence  limit  for  R  ,  then: 

s 


PLX  >  H 


s,L(a) 


]  =  1  -  a 


(6) 


holds. 

Prom  equation  (6),  R  should  be  the  (1-a)  percentile 
point  of  the  simulated  distribution  of  R_  T(a\. 


o—o 


For  the  investigation  a  will  be  chosen  as  .20  and  An    v 
will  be  defined  as  the  (l»a)  percentile  point  of  the  distri- 
bution  of  Rg  L ,    y      In  order  to  investigate  the  distribution 
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of  the  random  variable,  a  simulation  by  digital  computer  is 
used  to  generate  the  values  needed.   Five  hundred  values  of 

A 

Rg  T(a\  will  be  calculated  and  ordered  such  that: 

ks,Ua)^<K,Ua){2)<K,LM0)   •=  •  •  *   *s,L(a)(500> 
Since  the  value  of  a  chosen  for  the  simulation  is  .20,  the 

A 

value  of  R  _ ,    x  will  be  found  such  that  20$  of  the  ordered 
s,L(a)  r 

th 
values  are  above  it.   Thus  the  80   percentile  of  the  distri- 

A 

bution  of  R   r /„  \ ,  A/n    \,  is  found  and  should  equal  R_. 

The  measure  of  accuracy  will  therefore  be  the  value 

|r„  -  A/,    \l  .   The  accuracy  of  the  method  will  be  checked 
Is    (1  -  a)  | 

for  different  combinations  of  the  failure  rate  distributions 
among  the  four  components.   The  sample  size,  N*,  to  be  tested 
will  consist  of  50,  10,  and  500  units.   The  mission  times 
for  each  component,  T  ,  ,  will  be: 

(1)  .5  mission  units 

(2)  5.0  mission  units 

(3)  To.  such  that  -  J^m[B(T0l)]  =  Z(l)    (?) 

oi 

where  Z(l)  =  |^j  (8) 


The  third  test  was  selected  as  an  ad  hoc  rule  because  of  the 
intuitive  feeling  that  this  time  would  give  a  good  measure 

A 

of  accuracy  to  the  procedure.   The  distribution  of  Rg  T(a\ 
will  therefore  be  dependent  upon: 
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(1)  a 

(2)  H±        i  =  1,2,3,4 

(3)  T  A       i  =  1,2,3,4 

oi 

(4)  R±        i  =  1,2,3,4 

The  mean  and  the  standard  deviation  of  the  distribution 

A 

of  Rg  L(a)  will  also  be  calculated  in  order  to  better  eval- 
uate the  accuracy  of  the  model.   The  procedure  for  the 
computer  simulation  of  the  model  is: 


(1)  Generate  random  times,  T'..  i  =  1,2,3,4  j  =  l,2, N 

for  each  of  the  four  components  in  series. 

ml         ml         rp  I  m  I 

1  11    l   12    l   13  •  •••   •   -  ■   IN 

ml         ml         rp»  ml 

1  21    L   22    L   23  '   *   '    L   2N 

ml        mi        ml  ml 

1  31  L   32    33  '  '  '     3N 

mi        ml        mi  mi 

1  41    l   42   l   43  '   '   '     4n 


(2)  Test  these  times  against  the  selected  mission  time, 

T  . ,  i  =  1,2,3,4  and  if  T*4  .  <  T  x  ,  record  a  failure 
oi'  ij    oi' 

or  if  T'**  >  TQi ,  record  a  mission  success  and  then 
set 

in  i      -;  p  in     <  T 

1  ij   Xi  X  ij    ioi  I 


|   01 


V 


If  ."  ..  -■    [' 

1  J       'J  L 


(3)  From  the  equations  (1)  -  (6)  outlined  i  >   Chapter  II, 
calculate  the  resulting  values  of  Rg  rf„\   and-  Hs« 

(4)  Repeat  the  above  steps  five  hundred  times. 
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(5)  Order  the  computed  values  of  H  T,  x. 

(6)  Calculate  the  mean  and  standard  deviation  of  the 
resulting  distribution  and  find  the  80^  per- 
centile point,  An  _a\. 

The  complete  computer  simulation  is  found  in  Appendix  C.   The 
procedure  used  is  shown  above  and  comment  statements  have 
been  used  to  aid  in  following  the  steps  of  the  program. 


16 


CHAPTER  IV 

RESULTS 

For  the  procedure  stated  in  Chapter  III,  the  failure 
times  generated  will  be  either  gamma  distributed  or  trun- 
cated normal  distributed.   The  parameters  of  the  distri- 
butions were  chosen  so  that  the  reliability  of  the  component 
was  a  preselected  value  or  the  failure  rate  curves  would  be 
representative  of  those  curves  used  by  contractors.   The 
following  cases  were  simulated: 


CASE 

COMPONENT 

DISTRIBUTION 

PARAMETERS 

!i 

I 

1,2 

3,4 

gamma 
tr.  normal 

r=1.5,  s=.036 
/A=11.28,  <T=4 

.9950 

.9950 

II 

1.2 

3,4 

gamma 
tr.  normal 

r=1.5, 

s=.036 
<r=io 

.9950 
.9984 

III 

1,2 

3.4 

gamma 
tr.  normal 

r=1.5, 
m=-30, 

s=.036 
<r=45 

.9950 
.9707 

IV 

1,2, 

.3s 

,4 

gamma 

i*l. 5 t 

s=.036 

.9950 

V 

1,2, 

,3, 

,4 

tr.  normal 

A=8, 

<r=4o 

.9810 

VI 

1,2, 

.3, 

,4 

tr.  normal 

/iA=20 , 

<r=4o 

.9850 

VII 

1,2, 

,3, 

,4 

tr.  normal 

XA=-30, 

(T=45 

.9707 

Each  case  is  specified  by  the  type  of  distribution  and 
parameters  applied  to  the  four  components  in  series.   The 
results  in  Table  I  are  for  the  first  three  cases  using  all 
possible  combinations  of  sample  sizes  and  test  times.   The 
results  of  cases  IV,  V,  VI,  and  VII  are  given  in  Table  II. 
These  cases  are  different  from  the  first  three  in  that  mis- 
cellaneous test  times  and  sample  sizes  are  used  in  the  simu- 
lation. 
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The  first  five  columns  in  the  tables  are  self  explana- 
tory and  the  last  is  the  value  of  TT  which  is  represented  by 

the  formula  TT  =   2  N,  X.    Toi .   It  is  a  measure  of  the 

i  =  1  1      1 

amount  of  testing  required,  relative  to  the  component  un- 
reliability. 

Appendix  A  contains  discussions  of  the  gamma  and 
truncated  normal  failure  distributions.   Appendix  B  contains 
the  computation  of  the  ad  hoc  test  times  used  in  Table  I. 
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CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

A  partial  explanation  of  the  procedure  developed  in 
reference  1.  ,  which  is  used  to  compute  an  upper  limit  on 
failure  rate  estimates  and  lower  limits  on  reliability 
estimates,  has  been  presented  in  Chapter  II.   A  simulation 
of  this  procedure,  without  the  assumption  of  an  exponential 
failure  rate,  was  proposed  in  Chapter  III  and  the  results 
tabulated  in  Tables  I  and  II. 

Examination  of  the  results  of  the  first  three  cases 
shows  that  the  average  accuracy,  "D,  for  the  various  planned 
test  times,  PTT,  are  as  follows: 


PTT 

"D 

.5 

.0174 

5. 

.01^7 

ad  hoc 

.0047 

Many  more  simulations  would  have  to  be  conducted,  using 
various  parameters  of  the  failure  rate  distributions,  in 
order  to  establish  final  results .   But  from  the  cases 
presented  it  would  seem  that  the  test  time  selected  from  the 
ad  hoc  rule  stated  in  Chapter  III,  produced  the  best  accuracy 
regardless  of  sample  size.   From  Cases  IV  and  VII  the  average 
accuracy  using  the  ad  hoc  rule  for  test  times  was  .0085.   In 
Cases  V  and  VI  test  times  of  1.0,  5.0,  10.0,  15.0,  and  20.0 
mission  units  were  used  and  the  accuracy  varied  from  a  low  of 
.0000  to  a  high  of  .0071. 
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APPENDIX  A 
A.   Gamma  Distribution 


sr  tr-l  e-st 


f(t)  =  -a     rrc  —  (9) 

7(t)  -     f{t] 

z(t)  ~  "rTFT 

r  +.r  -  1   -  st 
2(t)  =  s  t     e  (10) 

_st  r~1  (st)1 

e  sc  Z 


i=  0 


l 


To  evaluate  the  parameters  and  meet  the  requirements  that 
H.  =  .995  the  following  procedure  was  used: 

R(t)  =  p[T>t] 

H.    =  fl(l) 

R.    =  P[T>1] 

R.    =  p[2sT>2s] 

Ri   =  P^2r>2s^ 
.995  -  PL^2r>2s3 

Therefore  for  a  given  value  of  the  parameter  r,  and  the  com- 
ponent reliability,  a  chi-square  table  look  up  would  give  the 
value  of  the  second  parameter  s,   For  r  =  1.5»  s  =  .  03&.   The 
specific  value  of  r  was  chosen  because  the  failure  rate  curve 
is  representative  of  those  curves  used  by  contractors.   A 
plot  of  equation  (10)  with  the  selected  parameters  is  siiown 
in  Figures  IA  through  IIIA. 
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B.   Truncated  Normal  Distribution 


i  ( t  -M  \ 2 

2   '  


f(t)  =  — ^r-  e      *"      b(|>(u)  (11) 


wnere     u  = 


<r 


b  - 


Z(t)  = 


f(t) 


(12) 


H(t) 


i  {t  ~M^2 
•      b   e^   ~^" 

Z(t)  =  <££l— -—- 5-  (13) 

i  ft  -M    ' 

(Tv'2tt  x- 


Z(u)  = 


0-  fr(u) 

1-  f  (u) 


For  Case  I,  H.  was  specified  at  „995  and  the  following  pro- 
cedure was  used  to  evaluate  the  parameters  <T  and  M   . 


R.  -  P(T>1) 

Hi  -  i("flT">"ar") 

R.  =  b(l-$(^4)) 
.995  =  b($(-  i-^-)) 
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In  Case  I,  (5" and  M, were  selected  to  be  k   and  11.28.   from 
equation  (12),  b  =  1.0024.   For  Gases  II  and  III  the  values 
of  the  parameters  were  chosen  such  that  tne  failure  rate 
curves  would  be  flatter  and  again  better  approximate  the 
failure  rate  curves  used  by  contractors,,   In  Case  III  the 
restriction  of  a  lower  component  reliability  was  also  placed 
on  the  parameters.   Parameters  for  Gases  V  and  VI  were 
selected  as  miscellaneous  values  in  order  to  give  more  com- 
plete results  to  the  analysis „   Plots  of  equation  (13)  with 
the  selected  parameters  are  shown  in  Figure  IA  through  IVA. 
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1* 
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FIGURE   IA 
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45678 
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Gamma  and  truncated  normal  failure  rate  curves  for 
Case  II. 

FIGURE  II  A. 


0-  N(x, -30,^5) 
l~|-  r(1.5,.036) 
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3   T~   5   T   7    8 

time  (mission  units) 

Gamma  and  truncated  normal  failure  rate  curves  for 
Case  III. 

Truncated  normal  failure  rate  curve  for  Case  VII. 

FIGURE  III  A 
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O  -   N(x,8,^0) 
Q-    N(x, 20,^0) 


1  2  3  Ur  5        T™        ?89 

time    (mission   units) 

Truncated  normal   failure   rate   curves    for   Cases   V 
and  VI o 

FIGURE   IV  A 
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APPENDIX  B 

COMPUTATION  OF  AD  HOC  TEST  TIMES 

The  test  times  for  Cases  I,  II,  and  III  were  selected  as 
follows : 

1.  .5  mission  units 

2.  5<>  mission  units 

3.  Tn,  such  that  -  Ji_  ln(R(T   ))  =  Z(l)     (7) 

oi 

The  computation  of  the  planned  test  times  under  the  ad  hoc 
rule  for  the  distributions  used  is  given  below. 

A.   Gamma  Distribution 

R(t)  =  P(T>t) 

substituting  into  equation  (7) 


=i-  ln(R(T  .  ))  =  -  —  ln(P(T>T  .  )) 
T         oi      >v  oi 

oi  oi 


-i-  ln(P(2sT>2sT      )) 
oi 


-  ^-ln(H(T01))  =  -  J-ln(P(J(2r>2STol)) 
oi  oi 


Therefore  the  following  equation  must  be  solved  for  T 

oi 


-   1   ln(P(Jf|r  >  2sT   ))  =  Z(l) 
oi 


3^ 


Solving  equation  (10)  with  the  parameters  r  =  1.5  and 
s=.036,  Z(l)  =  .0076.   Substituting  into  the  above  equation 

-  In  (P(J(2  >  .072  TQl))  =  .0076  TQi 


Solving  gives  T  .  =2.4  mission  units 


B.   Truncated  Normal  Distribution 

H(t)  =  P(T>t) 
substituting  into  equation  (8): 


-  ■_!_  ln(R(T  .  ) )  =  -  JL  in  (P(T  >  T,  ) ) 

m  01  rn  Ol 

xoi  oi 


.  _i_  m  (R(T.  ))  =  -  JL.   ln(b(l-{  (Tol'M)) 
xoi  ioi  y 

Therefore  the  following  equation  must  be  solved  for  T  .  : 

-  in  (b(l-"$(  0l"M)))   =   Z(l)  Toi  (14) 

Solving  equation  (13)  with  the  parameters  for  Case  I, 
Z(l)  =  .0037.   Substituting  into  the  above  equation: 


T_  -11.28 


-  In  (1.0024  (1  -Cp(  01  )))  =  .0037  T 


oi 


Therefore  for  the  truncated  normal  distribution  in  Case  I, 

T   =1.0  mission  units.   Solving  equation  (14)  with  the 
oi 
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parameters  for  Cases  II  and  III,  the  following  times  result: 


Case  II  T   =  1.75  mission  units 
oi 


Case  III  T   =  1*2  mission  units 
oi 


36 


APPENDIX  C 

The  gamma  and  truncated  normal  failure  times  were  com- 
puted using  a  Monte  Carlo  method.   A  pseudo  uniform  random 
number  generator  was  used  to  generate  a  normal  (0,1)  random 
variate,  V.   The  random  number  generator  and  the  compiling 
of  the  normal  variate  was  programmed  in  machine  language 
and  was  used  as  a  subroutine  function  called  SNORM.   The 
random  number  generator  used  is  called  URN  and  is  available 
at  the  N.P.G.S.  computer  facility  library. 

The  gamma  distributed  failure  time  is  compiled  from  a 


chi-square  value  with  2r  degrees  of  freedom  and  a  parameter, 

2s 
to  the  sum  of  2r  squared  normal  (0,1)  variates,  the  gamma 


s,  such  that  the  gamma  variate  =  AfL_.   Since  X?.r   ^s  eclual 


distributed  time  is  compiled  in  the  computer  simulation  by 

2r    2 
2   V 

i  =  1 
the  value,  * .   The  normal  distributed  failure  times 

£S 

are  truncated  at  zero  and  they  are  computed  from  the  gener- 
ated normal  variates.   They  are  set  equal  to  (S~V  +  /t\, 
where  (Pand^are  the  parameters  of  the  distribution. 

The  general  flow  of  the  simulation  is  to  test  each 
component  seperately  for  N^  replications.   The  estimate  of 


number  586,  which  represents  equation  (1).   After  each  corn- 


that  component  failure  rate,  X.»  is  given  by  statement 


ponent  has  been  tested  N.  times  and  it's  failure  rate 


estimated,  the  system  failure  rate  is  given  by  Ao.   The 
point  estimate  of  system  reliability  is  then  computed  in 
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statement  591,  which  represents  equation  (2).   The  number 
of  failures  and  the  standard  deviations  for  the  trials  of 
all  four  components  are  compiled  in  statements  587,  and  588. 
The  resulting  formulae  for  the  lower  confidence  limits  of 
the  failure  rates  are  given  in  statements  616  and  62^  de- 
pending on  the  number  of  failures,  as  explained  in  equations 
(^a)  and  (^b)„   The  whole  process  is  repeated  five  hundred 
times  and  the  distribution  ordered.   The  mean  and  standard 

A  A 

deviation  for  Rs  T(a\  and  fi  ,  the  reliability  point  estimate, 
are  also  calculated  in  the  computer  simulation. 
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DIMENSION  M3).BETA(31)  ,      R (4 ) , RLHAT < 500 ) , 
*PSHAT(500),RSBAR(3),SDEV<3).Sn4),SIGMA<4), 
*TT(4,3> ,XMEAN<4),RLBAR(3),RSDEV<3) 

DATA  BETA/1.507, 1.369,1.309, 1.2 72 f 1.246, 1.227 f 
♦1.2 12, 1.200,1. 190, 1.1 81, 1.1 73, 1.167, 1.1 61, 
*1. 156, 1.153, 1.1 46, 1.1 44, 1.1 39, 1.135, 1.1 32, 
*1. 129,1. 126, 1.124, 1.1 22, 1.119, 1.117, 1.115, 
♦1.113,1.111,1.110,1.0/ 
C 

C      FEAD  IN  SELECTED  MISSION  TIMES  TT. 
C      READ  IN  SAMPLE  SIZE  JD. 

C      READ  IN  PARAMETERS  EOR  GAMMA  D ISTR I BUTION ,R , S I . 
C      READ  IN  PARAMETERS  FOR  TR.  NORMAL  DIST.  SIGMA, XMEAN, 
C      BETA  CORRECTION  FACTOR  IS  READ  IN  AS  DATA 
C 

JD  =  50 

READ  3,  ( (TT( I, J) ,1=1,4) ,J=1,3) 
READ  3,  (R( I ), 1=1,4) 
READ  3, (SI< I), 1=1,4) 
READ  3, (SIGMA( I ), 1=1,4) 
READ  3. (XMEANf I ), 1=1,4) 
3   F0RMAT(4F10.3) 
C 

C      M  IS  INDEX  FOR  500  CALCULATIONS  OF  RSL  . 
C      I  IS  INDEX  FOR  4  COMPONENTS  IN  SERIES. 
C      J  IS  INDEX  FOP  NUMBER  IN  SAMPLE  SIZE. 
C      K  IS  INDEX  FOR  TEST  TIMES. 
C 

CO  9000  K=l,3 
RLDUM=0.C 
PSDUM=0.C 
C 

C      GENERATE  500  VALUES  OF  RSL. 
C 

PO  650  M=l,500 
TLAM=0.0 
CUM=C.C 
ChAT=0.0 
FTOTE=O.C 
DO  5  9C  1=1,4 
F  =  O.C 
SUM=C.O 
CO  5B5  J=1,JD 
C 

C      IF  COMPONENT  I  hAS  GAMMA  DISTRIBUTION,  GO  TO 
C      STATEMENT  NUMBER  30D  FOR  GENERATION  OF  TIMES. 
C 

C      IF  COMPONENT  I  HAS  TRUNCATED  NORMAL  DISTRIBUTION 
C      GO  TO  STATEMENT  400  FOR  GENERATION  OF  TIMES. 
C 

C      IN  THIS  PROGRAM  COMPONENTS  1  AND  2  HAVE  GAMMA 
C      FAILURE  DISTRIBUTIONS  AND  COMPONENTS  3  AND  4 
C      HAVE  TRUNCATED  NORMAL  FAILURE  DISTRIBUTIONS. 
C 

IF  (1-2)  300,300,400 
3C0  IA=2.*R(1) 
TOTAL=0.C 
DO  330  IB=1 , IA 
V=SN0PM(1 ) 
330  TOTAL=TOTAL+(V*V) 

TP  =  T0TAL/(2.*SI(1  ) ) 
GO  TO  500 
400  CONTINUE 
402  V=SN0RM(1) 

TEST=SIGMA(4)*V«-XMEAN(4) 
IFCTEST-O. 0)402, 402, 405 
405  TP=TEST 

500  IF(TP-TT(I,K) )510,550,550 
510  T=TP 

F=F*1.0 
GO  TO  56C 
550  T=TT(I,K) 
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c 
c 
c 
c 

c 
c 


56C 
585 

C>P6 

537 

S88 

5  90 
5  91 


613 
S14 

615 

616 

620 


624 
625 


^5 


SUM=S 
COM  I 
XM=JD 
hLAM  = 
TLAM  = 
FTOTE 
CHAT  = 
PUM=D 
f  ONTI 
PSHAT 
MOTE 
IF(NT 
NTDTE 
CONTI 
IFCFT 
TOUM  = 
ULAM  = 
no  TO 
CHAT  = 
TOUM  = 
FOUM  = 
TDUM* 
ULAM  = 
PLHAT 
P  L  DU  M 
RSDU^ 
CONTI 


UM*T 
NUF 

(F/S 

TLAM 
=  FTO 
CHAT 
UM  +  1 
NUE 
(M|  = 
=  FTO 
PTE- 
=  31 
NtlE 
CTE- 
RETA 
TDUM 

625 
CHAT 

BFT 
SORT 
TDUM 
(2.* 

(M)  = 

=  PLO 
=  RSO 
NUE 


UM)M2.*XN/(2.*XN+1.  )  ) 

♦  HL4^ 

Tfc+P 

♦HI.  AM/SUM 

•O/SUM 

FXP(-TLAM) 

TE+l.C 

31  )  6 14, 6 14, 61 3 

O.C)  615,615,620 

(NTOTFP-.S^? 
♦TOUM*DUM/4#0 

/TLA** 

A(NTCTE)*.842 

(4.»TLAM*TDUM*TnUM-t  CHAT  +  Tnijw*  TDJM* 

"CHAT* CHAT  I 

TLAM  +  TDUM*TOllM*  CHAT  ♦PDl«M|  /  2.0 

EXP(-IILAM) 

IJM*RLHAT(M| 

UM+«SHAT(M) 


COMPUTATITM  f:r  mean 
DISTRIBUTION  FOR  RSL 


AMD  STANDARD  DEVIATION  OF  THR 


69C 


7Ci 


70  3 


706 
705 


FSBAR(K )=RSDU*V5G0. 

RL8AR(K)=RLDUM/500. 

SLPUM=G.C 

SDUM=C.O 

CO    690     11=1, 50C 

SLDUM=SLOUM*(RLHAT(  I  I  ) -RL 3 AP ( K )  )* « 

SDUM=SDUM+(PSHAT( II)-RSBAR(K))**2 

SDEV(K)=S0RT(SLDUM/5GC.) 

FSOEV(K)=SORT(SDUM/50( . ) 


STATEMENTS  70C  TO  705  APE  USEn  FOR  ORDERING 
500  CALCULATED  VALHFS  OF  RSL  AND  SFLECTING 
c?0    TH  VALUF  OF  THE  DISTRIBUTION. 

DO  7C5  NX=1,101 

B=PLHAT(NX) 

NA=NX*i 

TO  7f6  N«=NA,5C: 

IF(R-RLHAT(NB) »7C3,7( 

TEMP=PLHAT(NB) 

PLHAT(N8>=B 

h=TEMP 

CONTINUE 

CONTINUE 

MK)=B 


THE 
THE 


6,7C6 


CUTPUT  STATEMENTS 
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